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Abstract 

We study the behavior of the soliton solutions of the equation 

where W' e is a suitable nonlinear term which is singular for e = 0. We 
use the "strong" nonlinearity to obtain results on existence, shape, 
stability and dynamics of the soliton . The main result of this paper 
(Theorem [T]) shows that for e — > the orbit of our soliton approaches 
the orbit of a classical particle in a potential V{x). 
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1 Introduction 

Roughly speaking a solitary wave is a solution of a field equation whose 
energy travels as a localized packet and which preserves this localization in 
time. A soliton is a solitary wave which exhibits some strong form of stability 
so that it has a particle-like behavior. In this paper we study the dynamics 
of solitons arising in the nonlinear Schroedinger equation (NSE): 



where (p(t, x) can be considered as a wave and \l/ e (t, x) is our soliton: a bump 
of energy concentrated in a ball centered at the point q = q e (t) with radius 
R e — > (for e — > 0). Considering this decomposition, the solutions to our 
equation can be thought as a combination of a wave and a particle, e occurs 
in the equation as a parameter. The main purpose of this paper is to show 
that for e and <£>(0, x) sufficiently small, our soliton behaves as a classical 
particle in a potential V(x). More exactly, we prove that the decomposition 

holds for all times and the bump follows a dynamics which approaches 
the dynamics of a pointwise particle moving under the action of the potential 
V (Theorem [T]); in particular the position q e (t) of the soliton approaches the 
position of the particle uniformly on bounded time intervals (Corollary [2]). 

The attention of the mathematical community on the dynamics of soliton 
of NSE began with the pioneering paper of Bronski and Jerrard then 
Frohlich, Gustafson, Jonsson, and Sigal faced this problem using a different 
approach [T5] . In the last years, several others works appeared following the 



first approach ([21], [25], [29], [30], [31]) or the second one ([I], [2], [3], [I], 

UM- & EE]). 




— + \w' e (V>) + V(x)ip. 
'm 2 



(1) 



A solution of our equation can be written as follows 



ip(t,x) = ^ £ (t,x) + (f(t,x) 



(2) 
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In this paper, we have studied equation (TjQ) which gives a different problem 
with respect to the ones mentioned above. Actually, in equation (0Q), unlikely 
the other papers, the parameter e appears in the nonlinear term and it will be 
chosen in such a way that \ip\ 2 approaches the delta-measure as e — > 0. Thus 
equation (pQ) describes the dynamics of a soliton when its support is small with 
respect to the other relevant elements (namely V(x), the initial conditions 
and its L 2 -size). Also the method employed here is different from those 
of the paper quoted above and it exploits and develops some ideas of [9]. 
Basically, we use the "strong" nonlinearity to obtain results on existence, 
shape, stability and dynamics of the soliton (Theorem [TJ). Finally, we notice 
that this method applies to a large class of nonlinearities and we do not make 
any assumption on the nondegenaracy or the uniqueness of the ground state 
solution (see the discussion in Section [TT^]) . 

1.1 Notations 

In the next we will use the following notations: 



Re(,2), lm(z) are the real and the imaginary part of z 



B p (x ) = B(x ,p) = {x G 



\x - x \ < p} 



B(x ,p) 



c 



N B(x ,p) 



S^iueH 1 : \\u\ 



L 2 



J{u) 



j [l |V«| 2 + W(uU dx 



J c = {ueH l : J (u) < c} 
d a V(x) 



doc ^ • • • doc <i 
J{u) = c 



where a£N, ii, . . . , i a G {1, . . . , iV} 



\d a V(x)\ = sup 

I a 2 = inf 

ueH 1 , J u 2 =a 2 

| • | is the euclidean norm both of a vector or of a matrix 

r = < U G H 1 , J(U) = inf J(u) > is the set of ground state solutions 
I 11^=1 J 



1.2 Statement of the problem 

First, we focus on the "concentration" properties of a soliton solution of eq. 
(pO) without the potential term V. We consider the following Cauchy problem 
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relative to the NSE: 

ip (0, x) =U e (x- q) e ip - x (4) 
where, with some abuse of notation, we have set 

WsW = J r2 W(e^M); W>W = J^W>(e^W)^- y (5) 

and W : M + — > R is a real function which satisfies suitable assumptions (see 
(ii) below). U denotes a ground state solution of the equation 

-.LaU + ]-W'(U)=u> 1 U (6) 
2m 2 



namely a function such that 

IWI 



2m 



J dx = c 



with 



IVtt 



2 



c = inf / J L + cfe (7) 



It is well known that we can choose U radially symmetric and positive. 
Direct computations show that, by virtue of (jSJ), the function 



—u( x - 

satisfies the equation 



- ^AU e + l -W' £ {U £ ) = u e U 6 . (8) 

where 

U e = -T- 

Moreover U E (x) is a ground state solution of (jHJ). In many cases, the ground 
state solution U is unique up to translations and change of sign, but we do 
not need this assumption. 

Notice that the choice of W E given by implies that 
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for every e > 0. 

Direct computation shows that the solution of 
following soliton 



with 



if) a . (t, x) = UAx-q- vt) e*^*- 4 "*) where v = — 

m 



lo = lo £ + -tod' 



is given by the 
(9) 



Thus if) £ (t, x) behaves as a particle of "radius" e living in the point 



Since ||V>,, E (*, 



q = q + v t 

| 2 = 1 for every e > 0, if e — > 0, we have that 
|^ g>£ (t,x)| 2 -> 5(x-q-vt) in V (R N ) Vt G R, 



(10) 



where 5 (x — x ) denotes the Dirac measure concentrated in the point Xq. 
The energy E £ (if>) of the configuration ip is given by 



so the energy of if) q E is 



2m 



dx. 



IV/7 I 2 \ I 



Thus if) qs (t, x) behaves as a particle of mass m: p can be interpreted as 
its momentum, ^p 2 = \mv 2 as its kinetic energy and 



+ W £ (U £ ) dx 



Co 



2m 

as the internal energy; here Cq is a constant defined as follows 

IVC/if 



c := 



2m 



+ W(£/i) da; 



The aim of this paper is to study the dynamics of the solitons in the 
presence of a potential V(x) namely to investigate the problem 



-±-Aif) + \W' e {iP) + V(x)i/j 



ip(0,x) = ipoJx) 



(P) 



where ip s satisfies the following assumptions 

*P 0iE (x) = U £ (x-q)e l ^ + ^ £ (x), cp , E e H\R N ); (12) 
= 1 (13) 

£ £ (Vv) < % + M ( 14 ) 

with M > independent of e; here _E e (-0) denotes the energy in the presence 
of the potential V : 



2m 



dx. 



1.3 The main results 

We make the following assumptions: 

(i) W : R + — > R is a C 2 function which satisfies the following assumptions: 

W(0) = W'(0) = W"(0) = (W ) 

\W"(s)\ < Cl |s| g ~ 2 + c 2 \s\ p - 2 for some 2 < q < p < 2* (W x ) 

4 

W(s) > -c\s\ u , c > 0, 2<u<2 + — for s large (W 2 ) 

3s G R + such that W(s ) < (W 3 ) 

(ii) V : R^ — >■ R is a C 2 -function with bounded derivatives which satisfies 
the following assumptions: 

< V(x) <V <oo; (V ) 

The main result of this paper is the following theorem which describes 
the shape and the dynamics of the soliton \l/ £ (t, x): 

Theorem 1 Assume (i) and (ii); then the solution of problem / fP]) has the 
following form 

i/> e (t,x) = * e (t,x) + <p e (t,x) (15) 

where ^ £ (t,x) is a function having support in a ball Bn E (q £ ), with radius 
R e — > and center q £ = q £ (t). Moreover, 

\\\V e (t,x)\-U £ (x-q £ (t))\\ L2 ^0 ase^O (16) 
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uniformly in t, where U £ = -^j^U (|) and U is a ground state solution of 
The dynamics is given by the following equations: 



with initial data 



p £ (t) = -W(<fe(t)) + F £ (q £ ) + H £ (t) 



q £ (0) = q + o(l) 
p £ (0)=p + o(l) 



(17) 



(18) 



where 



(a) q £ {t) is the bary center of the soliton and it has the following form: 



x \if?J dx 



|*e| dx 



(b) m £ (t) = m J \iff £ (t, x)\ 2 dx = m + o(l) can be interpreted as the 
mass of the soliton, 

(c) p £ {t) is the momentum of the soliton and it has the following 
form: 

p £ {t) = lm J (t, x) ^ £ (t,x) dx; 

(d) K £ (t) and H £ (t) are errors due to the fact that the soliton is not a 
point and 

sup (\H £ (t) \ + \K £ (t)\) ^ as e^O. 

(e) F £ (q £ ) is the force due to the pressure of the wave ip £ on the soliton 
and F £ — > in the space of distributions, more exactly we have that 



Vr ,ri, 



F £ {q £ ) dt 



<c(e)(l + \r 1 -r \) 



where c(e) — » as e — > 0. 
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Corollary 2 Let q and p be the solution of the following Cauchy problem: 

m = i Pit) 

(19) 

p(t) = -W(q(t)) 
with initial data 

q(0) = g e (0) , 9n , 
p(o)= Pe (o) [Z(j) 

where q £ (t) and p e (t) are as in Th. d Then, as e — > 

(q £ (t),p £ (t))^{q{t),p{t)) 

uniformly on compact sets. 

Let us discuss the set of our assumptions. 

Remark 3 The conditions HWp\) and jV^ ) are assumed for simplicity; in fact 
they can be weakened as follows 

W(0) = W'(0) = 0, W"(0) = E 

and 

Ex < V(x) < Voo < +oo. 

In fact, in the general case, the solution of the Schroedinger equation is mod- 
ified only by a phase factor. 

Remark 4 In FJ|/ the authors prove that if (ii) holds equation (T3|) admits 
orbitally stable solitary waves having the form Op. In particular the authors 
show that, under assumptions $Wq\ ), jWi\ ), and jW^D , for any a there 

exists a minimizer U(x) = U a (x) of the functional 

J{u)= [ \-\Vu\ 2 + W{u) ) dx 



on the manifold S a := {u G H 1 , \\u\\ L 2 = a}. Such a minimizer satisfies 
eq. (0>P where lo\ is a Lagrange multiplier. 

Remark 5 By our assumptions, the problem [F\l has a unique solution 

V> e C°(R, H 2 (R N )) n C\R, L 2 (R N )). (21) 
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Let us recall a result on the global existence of solutions of the Cauchy problem 
(EP (see /IS ESI [If j. Assume (TO , (TO and /or W . Let D(A) (resp. 
D{A 1 / 2 )) denote the domain of the self-adjoint operator A (resp. A 1 / 2 ) where 



A = —A + V : L (R J ) — )■ L , 

//K>0JeC 2 and |<9 2 1/| e L°° and tae imtta/ data ^(0,x) e D(A 1/2 ) 
£nen inere exists the global solution ip of IP]) and 

V>(*, x) e C° (R, D(A 1/2 )) n C^R, F-^R^)). 
Furthermore, ifip(0,x) e -D(-A) taen 

V»(t,a;) G C (M, J D(A))nC 1 (M,L 2 (M JV )). 

In tais case, since D(A) C //^(R^), ( [Hp is satisfied. 



2 The dynamics of solitons 

In this section, after having stated the main dynamical properties of our sys- 
tem (Subsection I2. II) . we define explicitely the splitting (j2J) and the equations 
which rules the dynamics of the soliton. In all this section, it is not necessary 
to assume e to be small. 



2.1 First integrals of NSE 



Equation (b?J) is the Euler-Lagrange equation relative to the Lagrangian den- 
sity 

,2 ,, , , , .!-> ^2) 



1 



C = Re(id t W) - — IWI - W e ty) - V(x) 
2m 

Sometimes it is useful to write ip in polar form 

if>{t,x) =u{t,x)e iS{t ' x) . 



(23) 



Thus the state of the system ip is uniquely defined by the couple of variables 
(n, S) G R + x [R/ (27rZ)]. Using these variables, the action S =j Cdxdt 
takes the form 



S(u, S) 



— |Vn|" + W e (u) + ( d t S + J- \VS\ 2 + V(x) ) u 



2m 



2m 



and equation ([Pj splits in the two following equations: 



dxdt 
(24) 
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- ^-Au + W £ (u) + (d t S + ^- \VS\ 2 + V{x)j u = (25) 

ft(u 2 )+V.(«^)=0 (26) 

Noether's theorem states that any invariance for a one-parameter group 
of the Lagrangian implies the existence of an integral of motion (see e.g. [T7] 
or [?]). They are derived by a continuity equation. 

Now we describe the first integrals which will be relevant for this paper, 
namely the energy, the "hylenic charge" and the momentum. 

Energy. The energy, by definition, is the quantity which is preserved by 
the time invariance of the Lagrangian; it has the following form 



A-\Vip\ 2 + w £ (ip) + v(x)\ip\ 2 



dx. (27) 



EM = 
Using (123]) we get: 

E e W) = f (J^ |V«| 2 + W £ (u)^jdx + J (J- \VS\ 2 + V(x)j u 2 dx (28) 

Thus the energy has two components: the internal energy (which, sometimes, 
is also called binding energy) 

J £ {u) = J (J^ \Vu\ 2 + W £ (u)^j dx (29) 

and the dynamical energy 

G(u,S) = j (^\VS\ 2 + V(x)^ju 2 dx (30) 

which is composed by the kinetic energy j- j \VS\ 2 v?dx and the potential 
energy J V(x)u 2 dx. By our assumptions, the internal energy is bounded 
from below and the dynamical energy is positive. 

Hylenic charge. Following [7] the hylenic charge, is defined as the 
quantity which is preserved by the invariance of the Lagrangian with respect 
to the action 

'ip i — y e %e ijj. 
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This invariance gives the continuity equation f l2T)j) . Thus, in this case, the 
charge is nothing else but the L 2 norm, namely: 

C(V>) = / \ip\ 2 dx = [ u 2 dx. 



Momentum. The momentum is constant in time if the Lagrangian is 
space-translation invariant; this happens when V is a constant. In general 
we have the equation 



d t (u 2 VS) = -u 2 VV + V • T 
where T is the stress tensor whose components are given by 

T jk = Re {d X] i)d Xk i)) 
-Sjk 



(31) 



(32) 



d Xj u d Xk u + (d Xj S d Xk S) u 2 - —S jk A (u 2 ) 



+5 jk [-W' e {u)u-W e {u) 



(33) 



Proof: It is well known that the stress tensor has the following form (see 
e.g. or 0) 

T jk = - Re [ - : , • d Xk i) ) + 5 jk C 



d(d x ^) 

Using equation (JP|) to eliminate id t ip, we get 



C = Re(id t ip4>) - — |Wf - W e (il>) - V(x) H 2 
2m 



Re 
Re 



l-\Vij\ 2 -W £ (ij)-V(x) 



-Aipijj 



1 ,_. l2 I- 



, 9/ JWf + 2^MV>-w 6 (v>). 



2m 

Moreover, we have that 

AW = A^i+A^ 2 = V-^iV^O+V-^V^-lVV'l' 



-A|^| 2 + -A|^ 2 | 2 - |W| 2 = -A|^| 2 - |VV| 4 
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Then 

c = -^AH 2 + \w' e (V) V> - wm 

We recall that, if z = a + ib is a complex number, by definition, ^ = 

da +l db- SO 

and the conclusion follows from direct computation. 
□ 

If V — const, eq. fl3Tl) is a continuity equation and the momentum 

P(^) = J U 2 VS dx = Jim tyVif) ^ (34) 

is constant in time. Notice that, by equation fT2B]) . the density of momentum 
n 2 VS* is nothing else than the flow of hylenic charge. 
Let us consider the soliton (jHJ); in this case we have 

where Co is defined by ([7]), 

= i 

and 

P(^)=P- 

Now, let us see the rescaling properties of the internal energy and the L 2 
norm of a function having the form 

: = ^ (f ) • 

We have 



Ml* = ~w v {-) dx = v (0 2 d£ = I MIL- 



and 
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\Vu\ 2 + WJu) 



dx 



2m 

J-|V U |S + ^W(^ M ) 



dx 



1 1 

1 1 

2m~e N + 2 



x 



fix 



Ar+2 



e 



W(v(0) 



e N d£ 



h l T^l^iOf + Wiv^d^^u). 



2.2 Definition of the soliton 

In this section we want describe a method to split a solution of eq. ([P]) in a 
wave and a soliton as in (j2J). 

If our solution has the following form 

1> e (t,x) = u e e iS < = [U £ (x-£(t))+w £ (t,x)]e iS ^ 

where w is sufficiently small, then a possible choice is to identify the soliton 
with U(x — ^{t))e lS ^' x ^ and the wave with w £ (t, x)e tS< ^' x \ However, we want 
to give a definition which localize the soliton, namely to assume the function 
ty £ (t, x) to have compact support in space. 

Roughly speaking, the soliton can be defined as the part of the field ip £ 
where some density function p e (t,x) is sufficiently large (e.g, after a suitable 
normalization, p £ (t,x) > 1). 

For the moment we do not define p £ (t, x) explicitely. We just require that 
p £ (t,x) satisfies the follwing assumptions: 



p £ e C l {R N+1 ) and p £ (t,x) ^ as \x\ 
p £ satisfies the continuity equation 

d t p £ + V • J Pe 



for some J Pe G C 1 (M 7V+1 ; 



oo 



(35) 



In order to fix the ideas you may think of p £ (t, smooth approxi- 

mation of u £ (t,x) 2 . An explicit definition of p £ (t,x) is given in Section f3T2l 
However, in other problems, it might be more useful to make different choices 
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of it such as the energy density. We have postponed the choice of p £ since 
the results of this section are independent of this choice. 
Next we set 

X £ (t,x) = vV iP e {t,x)) 

where 

( if s < 1 
(p ( s ) = } s - 1 if 1 < s < 2 
[ 1 if s > 2 

So we have that x £ {^i x ) = 1 where p £ (t,x) > 2 and x £ (t,x) = where 
Pe(^) x) < 1 : thus you may think of x e (i, as a sort of approximation of the 
characteristic function of the region occupied by the soliton. Finally, we set 

* e (t,x) = V> £ (M)x £ (36) 
x) = ^ e (*,a;)[l-x e ] (37) 

\l/ e (t,x) is the soliton and tp e (t,x) is the wave; the region 

£ M = {(t,x) Gl w+1 | Kp(t,rr) < 2} (38) 
= {(^)6R w+1 |0<x E (t,x)<l} (39) 

is the region where the soliton and the wave interact with each other; we will 
refer to it as the halo of the soliton. 



2.3 The equation of dynamics of the soliton 

Definition 6 We define the following quantities relative to the soliton: 

• the barycenter; 



Qe(t) = 



fx 




2 dx 




2 dx 



• the momentum: 

Pe {t) = J WS £ \^ £ \ 2 dx 

• the mass: 

m £ (t) — m J {^e^dx. 

Remark 7 Notice that the mass of the soliton m e (t) depends on t. The global 
mass is constant (namely m) but it is shared between the soliton and the wave 
whose mass is J u 2 £ [1 — x 2 e ] & x 
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The next theorem shows the relation between q £ (t) and p e (t) and their 
derivatives. 

Theorem 8 The following equations hold 

qe = + -i- / (x-qe) [u 2 e VS £ ■ Vp £ - V ■ J Pe ] dx (40) 

m e m e Js Eit 

Pe = - [w\V £ \ 2 dx- [ [T-Vp £ + u 2 £ VS £ (V-J^)]dx. (41) 

Remark 9 The term J s T ■ Vp £ dx represents the pressure of the wave on 
the soliton; if e — > and dY* £ j is sufficiently regular then 

I T ■ V ' p £ dx = l T ■ n da 

where a £ = {x \ p £ {x) = 1} and n is its outer normal. 

Proof of Th. [SJ We calculate the first derivative of the barycenter. 



We have 
and 



jx 




2 dx 




2 dx 



ifx\v £ \ 2 dx (fx\y £ \ 2 dx)dfm 2 dx) 

Jl^efdx (J\V £ \ 2 dx) 2 
&Jx\* e \ 2 dz iJlVefdx j(x-q £ (t))j t (\y £ \ 2 )dx 

fl^efdx qe{) fl^efdx " Il^efdx 

V X 2 = Vy? (p £ (t, x)) = up' (p £ (t, x)) Vp £ = I Se>< Vp £ 



J t X 2 = -£<P(P e (t,x)) = <p'(p £ (t,x))d t p e 



where Is £it is the characteristic function of E £jt . 
So, we have 
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. J(x-g £ (t))i( X 2 u 2 e )dx 

qe{ ) ~ im 2 dx 

J(x-q e (t))( X 2 iu 2 £ +u 2 £ f t x 2 ) dx 
J\^e\ 2 dx 

J r n(x - q £ (t))x 2 f t u 2 dx - f Sst (x - g c (t))V • J P dx 
Irn l^sfdx 

For the first term we use the continuity equation (126]) . We have 



(42) 



J{x- q £ (t))x 2 ^_u 2 dx = J{x- q £ {t))V ■ 



2 vs £ 

5 m 



- / ( M 2 VS £ ) X 2 ^ + - f(x-q £ (t))u 2 £ VS £ -Vx 2 dx 
m J v 7 m J 

/M /) + - / (x - q £ (t))u 2 VS £ ■ Vp £ dx. 



Concluding, we get the first equation of motion: 

. , . pJt) f Est ( x ~ &(*)) [u 2 e VS e • Vp e - V • J P J 

= + — 1 

m £ m £ 
Next, we will get the second one. We have that 

Pe = J (i^ S ^ X 2 dx + J u 2 VS £ ^ X 2 dx. (43) 

Now, using ( IBTj) we have that 

J (^-u 2 £ Vs)j X 2 dx = - J W |^ e | 2 dx + j V • T x 2 cfe 

= - J VV\^> £ \ 2 dx- J T- Vx 2 dx 

= - [ W\^ £ \ 2 dx- [ T-Vp £ dx. 



e.t 



The second term of eq. (1431) takes the form: 

! u 2 VS £ -^x 2 dx = - j u 2 VS £ (V • J Pe ) dx. 



16 



□ 



It is possible to give a "pictorial" interpretation to equations (jlDj) and 
(jH]l . We may assume that u 2 £ represents the density of a fluid; so the soliton 
is a bump of fluid particles which stick together and the halo S £jt can be 
regarded as the interface where the soliton and the wave might exchange 
particles, momentum and energy. 

Hence, 

• m E (t) is the mass of the soliton 

• is the velocity of the fluid particles and VS" £ is their momentum 

So each term of the equations ()4Up and (j4ip have the following interpre- 
tation 

• ^^p- is the average velocity of each particle; in fact 

Pe (t) _ fVS £ \V £ \ 2 dx _ j^ |fr £ | 2 dx 
m £ m £ J\^ £ \ 2 dx 

• the "halo term" ^- J s (x—q £ ) [u 2 VS £ ■ Vp e — V • J P J dx describes the 
change of the average velocity of the soliton due to the exchange of fluid 
particles 

• the term — J VV \ ^ £ \ 2 dx describes the volume force acting on the soli- 
ton 

• the term — J s T • Vp £ dx describes the surface force exerted by the 
wave on the soliton 

• the term — J E t u 2 £ V S £ (V • J Pe ) dx describes the change of the momen- 
tum of the soliton due to the exchange of fluid particles with the wave. 

3 The limit dynamics 

In this section, we analyze the dynamics of the soliton as e — > and we end 
proving the main theorem i.e. Th. [TJ 
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3.1 Analysis of the concentration point of the soliton 

If i[) e (t,x) is a solution of the problem ((P|), we say that q £ {t) is the concen- 
tration point of ip £ (t,x) if it minimizes the following quantity 

f(q) = \\W E (t,x)\-U £ (x-q)\\ 2 L2 . (44) 

It is easy to see that /(g) has a minimizer; of course, it might happen 
that it is not unique; in this case we denote by q £ (t) one of the minimizers of 
/ at the time t. 

Basically q £ (t) is a good candidate for the position of our soliton, but it 
cannot satisfy an equation of type (TTT1) since in general it is not uniquely 
defined and a fortiori is not differentiate. q £ {t) could be uniquely defined if 
we make assumptions on the non degeneracy of the ground state, but we do 
not like to make such assumptions since they are very hard to be verified and 
in general they do not hold. Actually the position of the soliton is supposed 
to be q e (t) as in Def. El However, as we will see, q £ is useful to recover some 
estimates on q £ . So, in this subsection we will analyze some properties of 
q £ (t). We start with a variant of a result contained in [6]. 

Lemma 10 Given u G H 1 , we define (if it exists) q G M. N to be a minimizer 
of the function 

9^ \\U{x-q)-u{x)\\\ 2 . 

For any r\ there exists a 8{rj) such that, if u G J Co+5 t r i) p g 1 ( see section llJ]) . 
q exists and it holds 

\\U(x-q)-u\\ m <ri (45) 
u 2 dx < rj (46) 

M. N ^B(q,R v ) 

where R v = —C\og(rj) and U G T. 

Proof: The proof of (1451) can be found in [6]. If U G Y, again by [B] we 
know that, for R sufficiently large, 

f U 2 (x)dx< [ CVe-^K 

J\x\>R J\x\>R 

Thus 

/ U 2 (x)dx = C 3 p^e-^dp = C 4 R N e- c * R < e~ c * R 

J\x\>R JR 
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where the Cj's are suitable positive constants. We remark that R does not 
depend on U. 

Now, it is sufficient to take R v > — ^ log(?y) and by (H5l) we obtain fj46l) . 
□ 

We define the set of admissible initial data as follows: 
B £ , M = fax) = U £ (x- g ) e ip °- x + <p(x) : E £ (V) < ^ + M and ||^|| L2 = l} 

Lemma 11 For even/ 77 > 0, there exists e = e(j]) > such that 

\4> £ {t,x)\ 2 dx<r) (47) 



R N \B(q e ,eRr,) 

where ip £ (t,x) is a solution of problem FF\) , with initial data in B e M and q £ 
is the concentration point ofi/j £ . 

Proof. By the conservation law, the energy E £ (ip £ (t,x)) is constant with 
respect to t. Then we have, by hypothesis on the initial datum 

E £ (4> £ (t,x)) = E £ (i; £ (0,x))<^ + M. 

Thus 

UM,x)) = E £ (iP(t,x))-G(iP(t,x)) 



\VS s (t,x)\ 
2m 



2 



+ V{x) 



u e (t,x) dx 



< E £ (^ £ (t,x))<^ + M (48) 

because V > 0. By rescaling the inequality (j4"8~l) . and setting y = x/e we get 

J(\e N ^ e {t,ey)\)<c + e 2 M (49) 

We choose e small such that e 2 M < 6(rj). Then e N/2 ip £ (t, ey) G J c o+ 5 M n S u 
and so applying Lemma [TUJ 

e"|V> e (*,ey)| 2 dy<»7 (50) 

Now, making the change of variable x = ey, we obtain the desired result. 
□ 
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Lemma 12 Ifip £ (t, x) is a solution of problem FF\) . with initial data in B £<M 
and e sufficiently small, then 

\if) e (t,x)\ 2 dx = r](e) (51) 

where T](e) — > as e — > 

Proof. First we prove that for every rj > 0, there exists S\ (rj) > such 
that, if ip £ (0,x) G B £l ^^ M , we have 



/ 



s(y) j_ 

M ' R? 



we get 



x)\ 2 dx < rj. 

Arguing as in the proof of Lemma [TT], if £i (77) < min 

(H7|) . At this point, since S\ {rj) < -i- we have that E\ (77) R, q < ^ e\ (77). 
Now set 

e(rj) =minei (C) . 
»?<C 

Clearly, e (n) is a non-increasing function (which might be discontinuous) 
and e (77) — > as r] — > 0. Then it has a "pseudoinverse" function 77(e) namely 
a function which is the inverse in the monotonicity points, which is discon- 
tinuous where e (rj) is constant and constant where e (rj) is discontinuous. 
Moreover 77(e) as e — > 
□ 

3.2 Definition of the density p £ 

First of all we notice that, in Lemma fT2| it is not restrictive to assume that 

rj = r](e) > e. (52) 

Now we set 

p £ (t, x) = a(x) * u(t, x) 2 

where, a £ (s) G C°°, 



3 \s\ < 7/8 1 1 — 77 
|s| > 77^ (l + 77 



and 

\Va £ {s)\<r)-l (53) 
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Lemma 13 Take ip £ a solution of W\) with initial data in B eM . 

If \x — q e (t)\ <T}*(l — 2r}»j then p £ (t, x) > 3 (1 — rj) 
if \x — q £ {t)\ > 77s M+ 277s) then p £ (t, x) < 3r/ 

where rj = 77(e) as in LemmaVLB, In particular we have that 

£ M C B (g 6 (t),7j£ (l + 2^)) \ B (&(*), 77s (l " 2^)) 

where E £ji is defined by / fffgj) . 

Proof. If \x — q £ \ < 7/5 ^1 — 277s J , then 

|x — g e | + < fl — 2?7^J + ^fr\ < 77^ fl — 77^ 

and hence 

S(g e , y/e) C B (x, 77s (l - 77^ 
Then, by using Lemma [12], 

«,(«,,) = fa,J, y - x) u M >s[ u M 
J JB(x, n^/z-r) 1 / 4 ) 

> 3 / u £ (t,y) 2 dy> 3 (1 - 77) . 

If |x - g e (*)| > 771 (l + 277IJ , 

\x — q £ \ — \fe > 775 ^1 + 277M — ^fr\ > 77? ^1 + 77^ 

and so 

B (x, 77* (l + n*Y\ C R N \ v^). 
Then, using again Lemma [T2| 

J Jb(x, r)S(l+rfSjJ 

< 3 / u £ (t,y) 2 dy < 3i] 

□ 
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Clearly, p e = a £ * u\ G C 1 (IR 7V+1 ) and, by fl2"6"|) . it satisfies the continuity 
equation f l5o^) with 

J Pe = a £ * (u*VS e ) . (55) 

Therefore, the results of Section |2] hold. In particular, we have that the 
support of ty e (t,x) is contained in B (q £ , rj^ ^1 H- 277^^ ^ when r\ is suffi- 
ciently small (namely 77 < 1/3). Moreover, by (1541) . we see that the size of 
the halo is an infinitesimal of higher order with respect to the diameter of 
the soliton. 



3.3 The equation of dynamics as e — > 
Theorem 14 The following equations hold 

p, = -VV(q,(t))+F c (q c )+H c (t) 



where 



and 



sup (\H e (t)\ + \K e (t)\) as e^O 



FMe 



T ■ V ' p £ dx. 



Moreover we have that 
Vt ,ti, 

where c(e) — >■ as £ — > 
Proof. We set 



F £ {qe) dt 



TO 



^^(l + ln-rol) 



Keit) = — f (x-q £ ) [u 2 £ VS £ ■ V Pe - V • J Pe ] dx, 

Hx, £ (t) = [ u 2 VS £ (V • J Pe ) dx, 

H 2>£ (t)=VV(q £ (t))- J VV(x)\V £ \ 2 dx, 
H £ (t) = H h£ (t) + H 2<£ (t), 



(56) 
(57) 
(58) 

(59) 
(60) 
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and hence, by Th. El we need just to prove 

We estimate each individual term of K £ . We have that 

sup \x - q £ \ < 2 (r] 1/8 + 2r/ /4 ) < 3r] 1/s 

since q e (t),x G B(q e , r/ 1 ^ 8 + 2// 1 / 4 ). 

Also, by fl53l) and well known properties on convolutions, 

sup |Vp e | < sup \Va £ (x) * u £ (t, x) 2 \ 



< HVaJI roo • WuJ\ 2 T2 < 



If ^ 6 (0,aO G B e>w , by (J3DJ), we have 



so, by Lemma [121 



m 2 |V5 £ | < / u 2 £ 

< 772 [2mG(^)] 2 < const. rj l 



and in particular 

/ m 2 |VS £ | rfx < r/^ [2mG(^)] 5 < const.rj^ 



Finally, by (ESJ) 



sup | V ■ J Pe | < sup I (V ■ a e ) * (u^ViS's) 



< l|V-a e || Loo • / u 2 \VS £ \ 



< ,)? = const. 7}'*. 



By §1, 
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|S 6 , t | < |s(g e (t),7 7 i(l+27 ? i))|-|s(g e (t) J »7s(l-27;s 



r/s 1 + 2775 



A? 



Wat 



778 f 1 — 2/7 



A' 



JV 

u>nV 8 



1 + 2t]s) - ( 1 - 2t]b 



JV+1 



< unV 8 ■ 5Ntj 8 < const. n « . 
So, by <m,--§J 



\K £ (t)\ < [ \{x-q £ )[u 2 VS £ -Vp £ -V -J PE ]\dx 



< sup \x — q £ \ 



< sup I re — q £ \ 



u 2 £ VS £ ■ VpJ dx + 



[ |(v-Jj| 



sup |Vp e | • / |itgV5 e | + sup |V-J, 



2 

£ 1 / I "*e 
i6E E|1 JE E , t 



< 3?7 1/8 [const. n' 1 ^ ■ n 1/2 + const.n' 1/4 ■ |E e , t |] 

< const, r/ 1//s . 



< const.n 1 / 8 
Then, by Lemma [T2l 



^-1/4.^1/2^-1/4.^ 



\K £ (t)\^0 



uniformly in t. 

Now, let us estimate \Hi >e (t) \ ; by (|66i) and (|65|) we have 



\H h£ (t)\ < [ \u 2 £ VS £ (V-J PE )\dx 



< sup V • Jp 



\uivs E \ 



< const.— ■ T] 2 = const. rj 



1/4 



7/ 



By the above estimate, 
We recall that 



1 1* £ | 2 = l-o(l) 
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when e — y 0, and that supp\I/ £ C B(q e , r/ 1 / 8 + 2r] 1 ^). We have 
W{q £ (t)) = (1 + o(l)) y W( ?£ (t)) |* £ | 2 



and so 
\H 2 (t)\ 



= J VV(x)\y £ \ 2 dx-VV(q £ (t)) 

= J \W{x)- W(g £ )||vI/ £ | 2 cfe + o(l) j W{q e {t))\^ £ \ 

< \\V"\\c°<r») J \x-q £ \\^e\ 2 dx + o{\) J W(<fe(t))|* e 

< o(l) (||HIc°(r") + \\VV\\ c o {mN) ) = 0(1) 



for all t. By the above equation, (1681) and (1691) . the (J58J) follows. 

Let P = P (ip £ ) be defined by (1341) . By the definitions of p £ , and (JMj) . for 

every t G IR, we have that 



[P.(*)-P(*)l 



(70) 



and since T G L 



P. 



y v5(i^ £ i 2 -?i 2 ) 

|V£| u\ dx = cb(e) 

where ci(e) -» as £ -» 0. By (13TD 

P=y (-m 2 W + V • T) <£r 

/ u 2 £ W dx. So, by gTJ and (JSHD 
/ W (it 2 - |vl/ £ | 2 ) dx - / [T • Vp £ + « 2 V5 £ (V • J Pe )] 
= fw {u 2 - |^ £ | 2 ) dx + F £ (q £ ) - f u 2 e VS £ (V • J Pe ) dx. 
Then, by (IBTJl) and Lemma [Ti2"l 

/ u 2 VS £ (V • J P J dx- fw (u 2 £ - |* £ | 2 ) dx 



dx 



F E (q £ ) - (p e - P) 













< (i) + 1| vvii^. y |^-i* 

< o(l) + ||W|U / 

JR 
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B(q s , s/e) 



u 2 £ dx = Ci(e) 



where ci(e) — > as e — > 0. Finally by (|70]h Vr ,Ti, 



F e (q £ ) dt 



TO 



ri 



P £ -P dt 



+ ci(e) (ri - r ) 



< |p £ (n) - P(n)| + [p e (r ) - P(r )| + Cl (e) (n - r ) 

< 2c (e) + ci(e) (n - r ) < c(e) (1 + |r x - r |) 



with a suitable choice of c(e). 
□ 

Collecting the previous results, we get our main theorem and Cor. [2j 

Proof of Th. CD By the def. (i3B]l . (j57jl . Lemma [JJ and Th. EJ Theorem 
[TJ holds with 

□ 

Proof of Cor. [21 We rewrite ( 1T7|) . ( ITU]) and (}2"0]) in integral form and 
we get 



p e (*) = p 6 (0) - P Q VV(q £ (s))ds + f* [F £ (q £ ) + H £ (s)} ds 

q(*)=&(o)+ j;?^ 

p(t)= A (0)-jJW(q(*))cfa 
and hence, for any |£| < T 



(71) 
(72) 



\qe(t)-q(t)\ < 



Pe{S) 



m 



ds+ / \K e (s)\ds 



< U [ \p £ (s) -p(s)\ ds + a^e) 
Jo 

where, by (I55|) . ai(e) — > as e — > and 

\p £ (t)-p(t)\ < [ \VV(q £ (s))-VV(q(s))\ds + [ F £ (q £ )ds + [ \H £ (s)\ds 
Jo Jo Jo 

< L 2 \q e (s) -q(s)| ds + a 2 (e) 
Jo 

where, by (|55|) . a 2 (e) -)■ as £ — )• 0. Then, setting z e (t) = |g e (t) - q(i)| + 
|p e (i) — p(t) \ i we have 



z £ {t) < L j z £ (s) ds + a(e) 
26 



with a suitable choice of L and a(e). Now, by the Gronwall inequality, we 
have 

< a(e)e Lt 

and from here, we get the conclusion. 
□ 
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